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ABSTRACT 

An experimental and analytical investigation of heat transfer from a horizon- 
tal i thin, square plate inside of an enclosure has been carried out. Experimental 
results have been obtained from both the upward-facing and the downward-facing sides 
of the heated plate. Starting with the integrated momentum and energy equations, 
approximate solutions have been obtained for heat transfer in the laminar and the 
turbulent regime that correlate well with own experimental data and the results of 
other investigators. Radiative heat transfer correction was given a special atten- 
tion. Effects of the enclosure-related recirculation of the test fluid, as well 
as effects of simultaneous heat transfer on both sides of the plate, caused an early 
transition, and indicated a high level of internal turtulence. 


INTRODUCTION 


Heat transfer from horizontal, heated plates by far has not received the atten- 
tion accorded to other geometries. There appears still to exist a general scarcity 
of representative experimental data, in particular in the turbulent regime. As heat 
transfer by free convection, for moderate temperature differentials, is of the same 
order of magnitude as that by radiation, the need for radiation correction during 
experimentation cannot be disregarded. Therefore, the experimental results below 
have been carefully checked for the effects of radiation. Also, equations of the 
boundary- layer type are set up, describing free convection on both sides of a hori- 
zontal plate of finite dimensions, heated electrically. Similarity is assumed for 
both the velocity and temperature profiles. These assumptions are then justified 
by a close agreement between the results of the analysis and the experimental find- 
ings and the literature on the subject. The solutions obtained apply for laminar 
as well as for turbulent free convection. Due to highly sheltered conditions inside 
the chamber where the present experiments were carried out, undesirable disturbances 
like uncontrollable small air movements were eliminated. Ample space between the 
walls of the chamber and the test plate should have made the resulting free convec- 
tive flow essentially undisturbed by the presence of the confining walls. 


SYMBOLS 

A area of specimen, m 2 

1 

2 / 

a thermal diffusivxty, m /s 
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a constant in reference velocity equations, s 
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specific heat at constant pressure, j/(K kg) 

Grashof number, g3L^(T w - T^)/v^ , Gr^/l? , m ^ 
heat transfer coefficient, W/(m^K) 

thermal conductivity, W/(m K) ; empirical constant in expression for a 
numerical constant in expressions for the Rayleigh number 
pressure, N/m^ 

Prandtl number, v/a 

heat flux density, w/m ; heat flow rate, W 
coordinate, radial direction, m; radius of test plate, m 
gas constant, j/(kg K) 

Reynolds number, Ur/v 
Rayleigh number, Gr Pr 
thermodynamic temperature, K 

dr/dt, velocity ; reference velocity, r-direction, m/s 

O 

specific volume, nr/kg 

dz/dt, velocity ; reference velocity, z-direction (perp. to plate), m/s 
coordinate perpendicular to plate, m 
Nusselt number, hL/k 

—i 

cubic expansion coefficient, K 
boundary-layer thickness, m 
emissivity of specimen 
kinematic viscosity, m /s 
mass density, kg/m^ 

Stefan-Boltzmann constant, W/(K^ m^) 
reduced temperature, T - T , K 


In turbulent regime, u, U, w, and W represent mean-flow velocities 
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Subscripts! 

m condition where maximum velocity occurs in the polynomial expression 

w " at the wall of plate 

oo " far away from wall 

EXPERIMENTAL TEST SET-UP 


In order to make the intended free-convective experiments reasonably accurate, 
it was decided to test the specimens) that were to be used as the horizontal plates 
in the tests, independently for emissivity under high-vacuum conditions. Therefore, 
a high-vacuum chamber was built first, that could contain pressures as low as 1-0 4 
mm Hg, with the mean free path of molecules inside long enough to eliminate convec- 
tion, and to leave only radiation as the mode of heat transfer. Conduction along 
the specimen support wire, and the heater and the thermocouple wires was calculated 
to be less than 1 % of the total heater input. For free-convective measurements, at- 
mospheric pressure inside of chamber was reestablished. 


As the total area of specimens, A^ , represented only 2 . 1 % of the wall area of 
the vacuum chamber (Fig.l), the required emissivity could be calculated from the 
formula 4 u 

Q = € t a A,(T/ - t;) (l) 


The emissivities measured with this set-up are seen in Fig. 2 . They all appear to 
be of the right order of magnitude, when compared with results from the literature, 
and appear to be adequate for making radiative-loss corrections in free-convective 
measurements more realistic . 


BASIC EQUATIONS OF THE PROBLEM 

Let's now consider natural convection flow of a Newtonian, single species fluid 
over a horizontal surface in an axisymmetric flow, with the coordinates r and z. 

The flow is in the steady state, the velocity field consists of u = dr/dt and 
w = dz/dt, and the corresponding temperature field is T = T(r,z) . The acting sur- 
face force is the pressure, and the body force is the gravitational force, g, all 
acting on a horizontal plate 86.0 mm in diameter which, in turn, gives rise to con- 
vective flow considered positive in the direction towards the center of plate; at 
edge of plate, r = 0 . We assume that this flow is amenable to analytical treatment 
using the boundary-layer theory approach, together with Boussinesq approximation to- 
wards the treatment of the variable physical properties. Also, there is assumed si- 
milarity of velocity and temperature profiles. Then the continuity, momentum, and 
the energy equations assume the form 
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with the boundary conditions! 


at z=0, u=w-0; T = T - const; at the edge of plate , r = 0, and u = 0 ( 6 ) 

w 

at z = 6 , u = 0, T = T | du/dz = dT/dz = 0 (7) 


On differentiating eq . (4) with respect to r, and integrating it with respect to z, 
we have 

- - = gp f 5 r— (T - T ) dz ( 8 ) 

p dr J z dr' oo / ' ' 

This approach to free convection on a horizontal plate, apparently first outlined 
by Stewartson [l], for two-dimensional flow, is applied here to three-dimensional, 
axisymmetric flow. Thus, from the continuity equation, we have 


yield 


>Zfdu , u> . 

Equations ( 8 ) and (9) » on substitution in eq. (3)* after partial integration, 

( 10 ) 


pZ/CJU , U\ - 
" = - + da 


i i (r f o u2 dz) = 86 f o£ f z k (T ■ T J dz ] dz - v 


Equation (5) may be transformed in a similar fashion after a partial integration. 
Letting 0 = T - T , and assuming that 6 (hydrodynamic) « 6 (thermal), there results 


cpif(r / 6 uQ dz) = -k (|^) 
p r dr' J o ' 'dz / w 


Additionally, letting the velocity and temperature profiles be 


( 11 ) 


u = U (z/S ) (l - z/ 6 ) 2 ; 0 = 0 w (i - z/s ) 2 (l2) 

as proposed, for example, by Eckert and Drake for a related problem [ 2 ], and letting 

U = G 1 r m ; ^ = C 2 r 11 ( 13 ) 

a speedy solution of eqs. (l0) and (ll) is obtained. With 8 = l/T^, we get 

i ^ (rtJ 2 6 ) = 17.5 gg 96 d 6 /dr - 105 v(u/6), and (l4) 

i s (*»„«) = 60 “V 6 (15) 

which, after eq. (13) has been substituted, yield 
(2m + n + l)C 2 C 2 r 2 ^ 11 " 1 = 1 7-5 g 8 9^ C 2 r 2n_1 - 105 (C^Cg) r m_n v and ( 16 ) 

(m + n + l) C 1 C 2 r m+n-1 = 60 a r" n /C 2 ( 17 ) 


For similarity solutions, both sides of these equations must be independent of r; 
this is true if m = l/5 and n = 2/5 for the present geometry. Then eqs. (l 6 ) and 
(17) yield, with Gr* = Gr/L3 

C 1 = 2.978 a Pr 4 / 5 (9/l^+ Pr )" 2 /- 5 (Gr *) 2 / 5 and 
C 2 = 3,548 Pr_2//5 (9/ 14 + Pr ) 1 / 5 (Gr *) -1 / 5 


45 


from which, together with eg. (13) is calculated the local heat transfer coefficient 



(18) 


while the average heat transfer coefficient becomes 

hjj = O.705 k Pr 2 /- 5 (0.643 + Pr)” 1 /- 5 ^*) 1 /-^” 2 /- 5 (19) 

In order to obtain a formula for the average heat transfer coefficient 
based on the sidf L of an equivalent square plate, we let L 2 = rrR 2 , such that L 
(equivalent) = rr 2 R, hp = and Nup = hpL/k, or 

Nu l = 0.886 Pr 2 /- 5 (0.643 + Pr)” 1 /- 5 Gr^/- 5 (20) 

For air (Pr = O.72), this generates the Rayleigh number, Ra^= Gr^Pr, based on L, or 

Nu l = 0.780 Ra^ 5 (21) 


for a horizontal square plate heated on the top side only. With the present choice 
of coordinates, 6 = 0 at the edge of plate. In the case of a horizontal plate, 
heated at both sides, 6 / 0 at the edges. According to the literature on this sub- 
ject (cf. reference 3, for example), a horizontal plate heated on both sides has 
its heat transfer on the upper side appreciably increased by formation of a warm air 
"bubble" on its lower side. It appears in the present case that increasing the area 
of plate by about one fifth generates an "effective area" that shows realistically 
the effect of heating of the bottom side of plate on heat transfer on the upper side 
in the form of an "effective radius" , R e ff 

B 8ff ' 1 - 095 "actual < 22) 

to be used as the upper limit in the evaluation of the average heat transfer coeffi- 
cient 

h _ __ = 0.816 k Pr 2 /- 5 (0.643 + Pr)” 1 /- 5 Gr 1 /- 5 ^ (23) 

II f ©I I 

which, for example, for Pr = 0.72, results in 

Nu l = 0.902 Ra^/- 5 (24) 

The velocity and temperature profiles applicable here as per eqs.(6) and (7) 
satisfy the physical conditions of the problem at hand, eqs.(2) - (5)? and reduce to 
zero for z > 6 . They are also for z < 6 qualitatively what can be reasonably expec- 
ted in free convective flow under conditions of conservation of mass and energy, 
regardless of orientation of flow, with exception of a small area near the plate’s 
center, where an upward-pointing jet is generated, as has been already observed by 
Stewartson [l]. This effect will be accounted for below in a separate discussion. 
Other velocity and temperature profiles have been considered. It may be seen from 
the comparison with the experimental results (cf . discussion below) that eqs. (2l) 
and (24) come already reasonably close to physical reality. This is perhaps the 
best justification of the methods used above, and for the particular u and 0 used. 


Experience teaches us that in free convection inside an enclosure (cf . McAdams 
[5] , p. 182) turbulence may develop already at relatively low Grashof numbers. For 
the case of heated, horizontal plates, this is explainable by the destabilizing ef- 
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feet of the upward-directed buoyant force on the boundary- like layer formed by the 
fluid flowing in the direction tangential to the plate. For developed turbulent 
flow, the last term in eq. (lO) may be replaced by the semi-empirical expression 

0.0228 U 2 (v/U 6 ) P ( 6 / 6 m ) P J 2 “ P (25) 

and the right-hand side of eq. (ll) may be defined as 

0.0228 c p U 0 (v/U 6 ) P Pr _ 2 / 3 ( 6/6 ) P J 1_P ( 26 ) 

p w m ' ' 

while the velocity and temperature distribution may be written (cf . Eckert and Jack- 

son [4], and reference 2 , p. 32b) as 

u = UJ (z/ 6 ) 1 / q (l - z/ 6 ) s ; 9 = 0^1 - (z/d) 1 ^] (2?) 

with the exponents q, s, and t determined later. 

Since eqs. (25) and (26) were originally based on the maximum velocity in the 
boundary layer, occurring here at z = 6 , a normalizing scale factor is introduced 
in eq. (27), assuring that u/ U = 1 .0 when z = 6 m applies. It may be shown that u 
peaks out when z = 6 m = 6 /(q*s + l), while the nominal length 6 itself signifies to 
what extent the effect of heating from the plate penetrated the fluid above. The 
numerical constant in eq. ( 25 ) may be approximated by the formula 

C(p,q) = 1.013 (0.92 q + 2.30 ) _P<1 ( 28 ) 

based on results of Wieghardt (cf. Schlichting [ 6 ], p. 60l). With p = 4 , eq. (25) 
becomes the Blasius formula, linked to q = 7 in eqs. (27) and ( 28)5 when p = l/ 5 , 
q = 9 applies. It is felt that with these changes, eqs. (25) to ( 28 ) will be repre- 
sentative of turtulent flow in general, regardless of orientation. Furthermore, 
eqs. ( 27 ) and ( 28 ) satisfy the boundary conditions for u and 0 at z = 0 and z = 6 

in a satisfactory manner; with these substitutions, with the required integrations 

carried out, eqs. (l 0 ) and (ll) yield the relations 

J 2 I^l/r^ruA) = I 2 gg 0 w 6d6/dr - K 3 U 2 (v/U 6 ) P J 2_P and (29) 

J I 3 (l/r) f^(iU 6 ) = K 5 U(v/tJ 6 ) P Pr ' 2//3 J 1_P ( 30 ) 

Assuming again validity of eq. ( 13 ), eqs. (29) and ( 30 ) yield after a few 
transformations the expression for the local Nusselt number, for m = l/(2+3p), and 

2m ’ Nu r = JK 5 (K 2 /K 1 ) m (K 5 /K 4 ) m(l " 3p) Pr" m(l ' 3p) (l + Pr 2 ^/^ ) Aa r " ( 31 ) 

and the average Nusselt number for an equivalent square plate with side L, 

Nu l = JK 5 (K 2 /K 1 ) m (K 5 /K 4 ) m(l " 3p) Pr" m(l_3p) (l + Pr 2 ^/^ ) _m “ 3p) Ra L m 

(31a) 

In eqs. (3l) and (31a), the use of J should be considered entirely empirical; the 
constants below stand for the integrals based on eqs. (27) and ( 28 ): 

K 3 = K 5 = C(p,q) ( 6 / 6 m ) P 

= (2m + n + iy(U^6)j^u 2 dz= (2m+n+l) I 
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K 2 = n /(0/ d S/dx) • J^(J^ |^0 dz) dz = n Ig = n/(2t + l) 

= (m + n + l)/(U0 w ) * J^uQ dz = (m + n + l) 1^ 

With the use of the Blasius formula for wall shear stress and heat transfer, 
and s = 8 and t = 4 in eq. (27) results have been obtained that for 5x10-5 < Ra < 
10° deviate for less than ±1.5% from the least-squares fit curve representing the 
present experimental data, or with m = 4/ll « 0.364, and Pr = 0.72, 

Nu L = 0.104 Ra L 0, 364 (32) 

Here, m = 0.364 is somewhat lower than the experimental value of m = 0.384. The ex- 
perimental value of m was closely approximated with m= O.385, when p = l/5» <1=9. 
s = 8 and t = 1 in eqs. (27) and (28) was used, but the constant in eq. (31a) was 
6% too low vis-a-vis experimental data. In order to get the commonly cited in the 
literature value of m = l/j for the turbulent, free-convective heat transfer from 
the upward-pointing, horizontal plate by the present method, p = l/3» q = 5> s = 3» 
and t = 5 would have been required in eqs. (25) to (28), to yield 

Nu l = 0.150 Ra^ 3 (32a) 

Equation (32a) seems to fit experimental data of others, for free convection with- 
out recirculation, reasonably well. 

The analysis outlined above has been carried out along the lines suggested in 
references 2 and 4, but modified in several respects to fit the physical require- 
ments of a horizontal plate of finite dimensions, heated and pointing upward. As is 
stated in reference 4, the exponent of Ra, m = 2/5f fitted very well experimental 
data for turbulent heat transfer, for free convection on vertical plates and a ver- 
tical cylinder. It is seen that semi-empirical approaches can be used to advantage 
where the exact analysis still fails to give answers of practical significance. 


EFFECTS CF FREE CONVECTION ON HEAT TRANSFER AT CENTER OF A HORIZONTAL PLATE 


It is obvious from the consideration of the continuity equation alone that 
free convective flow, for the present geometry, will generate near the center of 
plate a vertical jet, known as the thermal plume. For purposes of the analysis be- 
low, we assume that temperature distribution within the plume is still given by 
eq. (12), but velocity follows the profile u = U(z/6 m ), z < 6 m , u = U, z > 6 m , and 

U = -a r j W = 2 1 a z (r = z = 0 at center of plate) (33) 


where i = 1, axisymmetric flow, and i = 0, plane flow. U and W in eq. (33) satisfy 
identically the continuity equation. By considering a volume element (2it r) 1 dr6 
and the heat fluxes in and out by convection, as well as by conduction at the bot- 
tom through the area element (2nr) 1 dr, after the substitutions for u and 6 have been 
made, there results the expression 


'V P d 


dr J o 


ue to) = k(||) M 


( 3*0 
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Taking 0 from eq. (12) and U as stated above, the integral in eq. (3*0 may be writ- 
ten, for 6 m /6 « 1 , as I = -ar^ 6 / 3 . It yields for 6 (which is actually the height 
of the present thermal plumej, the expression 

6 = ( 3 a/a ) 2 (33) 

Using eq. (l 8 ), we may write Nu r = 2R/6 . The still required expression for a 
in eq. (33) is now to be derived. If differential pressure at the center of the plate 
is proportional to k*+gp(30 W z/R , working through a distance dz, that accelerates a 
fluid element from velocity 0 to ¥ (here k_^ is the constant of proportionality) , one 


gets p k — £ 3g0 z dz = p dA j w (dw/dz) dz, on assumption of inviscid flow 

no 


W 


there, the simple result 

v = (ege i^/r ) 2 z 


From comparison with eq. (33) » eq.* (36) yields 

a = (pge^k^)^ = (k^ Grv 2 /at T / 2 , and 
A. LA. 

NUp = ( 2 / 3) 2 k Pr 2 Gr R * 

i. 

For air with Pr = 0.72, and L = tt 2 R, there results 


(36) 

(37) 

(38) 


Nu l = 0 .864 k Ra L * 


(39) 


Integral methods generate useful results that still depend on the approximating po- 
lynomials needed to represent u and 0 . Thus, if in the present case we use eq. (l2) 
for both u and 0 , tut with the reasonable constraint that, at z = 6 m , u = U, so that 
eq. (33) is satisfied, we obtain from eqs. ( 3*0 and ( 37 ) 

JL 

Nu l = 0.710 k Ra L 4 ( 39 a) 

Application of the same procedure to the case of an infinite strip (i = 0 in 
eqs. ( 33 ) and ( 3*0 • u and 0 according to eq, (12), and scale factor J = 27/*+ applied 
to u) , yields for a strip of the width 2 b 

Nu b = 0.62 js Ra b * (40) 

It is to be emphasized that eq. ( 3*0 > obtained from a direct summing up of the 
convective and conductive contributions to the control volume, differs in sign from 
eq. (ll), obtained from the exact integration of the energy equation, as the chosen 
coordinate system is different. Also, as the hydrodynamic information on 6 is not 
required, the momentum equation is not at all involved. This greatly simplifies the 
problem. The use of scale factor J here is needed to satisfy the continuity, eqs. (33)* 


RELATIONS FOR A HORIZONTAL DOWNWARD-FACING PLATE 

It appears that for this geometry the conventional boundary- layer approach is 
not possible (cf. reference ( 3 ), for example). In the present case, the only reason 
for the movement of fluid is its thermal expansion, due to the contact with the wall 
at T w . Therefore, the resulting flow is similar to that occurring in forced, axisym- 
metric stagnation flow, for which an exact solution of the Navier-Stokes equations 
exists (reference 6 , p. 100) . One can expect the fluid temperature near r = 0 
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(4i) 


to vary, on the average, at some distance below the plate's center, as 

0 = ( 0 ^ 2 ) (1 - r 2 /H 2 ) 


while simultaneously the velocity terms for the potential flow related to the pre- 
sent problem follow the relations satisfying identically the continuity equation: 

U = ar ; ¥ = -2^az (42) 


that applies to any incompressible flow, to any geometry. If the fluid in question 
is perfect gas, the term SP/dr in eq. (2) may be now transformed as ( with Pv = RT) 


3P _ /dP\ /3T\ _ p 
dr v 3r ® 


9 M r/E‘ 


or - 


3P 

dr 


0 r 
w 


using eq. (4l) for T = 0 + T . In quiescent gas, P = g/v times an element of length 
assumed here to be proportional to R . As a is still a free constant , we may assume 
to be that constant of proportionality, without the loss of generality. Then, 


- ~ gg0 r/(Rv) = k^ Gr v 2 3 

3r — w — 


<Bv) 


-1 ;letting a 2 = k^ Gr*v 2 /R, 


there results 


- H = P a 2 r (43) 

With the present definition of a, the pressure-gradient related term in the 
momentum equation is linked to the acting temperature differentia^, 0 W . Writing 
u = r f (z) , and w = -2 f(z), and using a new variable = ( a/v) J z , eq. (3) above is 
readily solved. This solution is considered as good for the complete Navier-Stokes 
equation (reference 6, p. 100) . The corresponding solution of the energy equation 
is available (cf. references 7 and 8). For Pr near unity, we have 


Nu r = O.763 Pr 0,4 Re r 2 (44) 

This is a boundary-layer type solution, based on a exact solution of the Navier- 
Stokes equations. As here Re r = Ur/v, wi^h a from eq. (43), we get Re r = kf • 

(Gr r r/R)l/2 , with Pr = O.72, and L = tt 2R , there results for air 

Nu l = O. 63O k Ra L * (45) 

where k is an experimental constant with a value k « 0.5, in view of the constant 
K = 0.27 and m = l/4 in the original correlation due to Saunders and Fishenden, for 
example (]9]» It must be stressed that the present approach is valid only in the re- 
gion where eqs. (4l) and (42) continue to apply, that is, near the center of a cir- 
cular plate (i = l), and near the center of a long strip, where r is the distance 
away from that center, (i = 0) . It is of some interest to note that a solution can 
also be obtained through a direct integration of eq. (ll) here, for the same veloci- 
ty and temperature distributions as in eq. (l2), but using a normalizing factor J = 
27/4, such that eq. (42) is satisfied. Then from 6 = (a/a)!/ 2, 2. 11, and Pr = 0.72 
we have the formula 

_i 

Nu^ = 1 .00 kj Ra^* (46) 


The term k' expresses the effect of geometry at hand as well as that of the approx- 
imating polynomials used, k' = 0.25 yields K that comes close to the above Saunders 
and Fishenden result. 

It is seen from the discussion above that methods, based both on the exact so- 
lution of the Navier-Stokes equations, and on the integrated energy equation method, 
lead to results that are experimentally verifiable . The relations for streamline 
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and temperature distribution, calculated by Miyamoto et al. flO] agree qualitatively 
with our eqs. (4l) and (42) and are also suggested by the scheme shown in Fig. 3 • 

In reference 10 are also reported results of Sugawara and Michiochi [11], where K = 
0.264 and m = l/4, that are based on the boundary- layer approximation. 

In addition to our eqs. (45) and (46) , and to solutions reported in refer- 
ences 10 and 11, there exist also several approximate analytical solutions for the 
present geometry with m = l/5 proposed (cf. references 12 and 13)- As an expla- 
nation for this one should consider that as the exact solution of the governing mo- 
mentum equation is only valid near r = 0, and with m = l/4 applicable for the associ- 
ated heat transfer problem, the flow starting far away from the stagnation point, 
at the edge of the plate, would result in a quasi-boundary layer flow for which, at 
least initially, m = l/5 would be entirely appropriate. The real flow must be actu- 
ally nonsimilar for a good portion of 0 < r < R. 


DISCUSSION OF EXPERIMENTAL RESULTS 


It is interesting to note that the present experimental results show that 
transition to turbulent flow has already occured as early as at Gr = 5x1 0?, in a- 
greement with the discussion by McAdams [5] » P-182, as the side effect of the con- 
finement in the relatively small vacuum chamber. Our own experimental results on 
free convection on an upward-facing, horizontal heated plate are shown in Fig. 4. 

It is seen from Fig. 4 that, for Ra< 6 x 10^ , the experimental data show a rea- 
sonable agreement with eq. (24) here, but the number of experimental points for 
what is considered laminar free convection on a typical upward-pointing, horizontal 
plate, heated on both sides, is not sufficient for a definitive finding. For com- 
parison, let's consider the formula proposed by Lewandowski and Kubski [l4] for a 
horizontal plate 100 x 60 mm, heated on one side only, 

Nu = 0.66 Ra 1 /- 5 (47) 

4 7 *t 

valid for 10 < Ra < 10 , with a ± 10% error margin, based on experiments with dis- 

tilled water, glycerine, and soybean oil. Equation (47) is very close to our 
eq. (2l) above. On the other hand, Yousef et al. [l5] , using Mach-Zehnder inter- 
ferometer, and for lCr<Ra <10?, in air, obtained 

Nu = 0.622 Ra 1//4 (48) 

for square plates with L = 100, 200, and 400 mm. A similar exponent of Ra has been 
also obtained by a number of other investigators, for air and for the same geometry: 
Fishenden and Saunders [l6] , Bosworth [l 7] , Al Arabi and El-Riedi [l 8] , for ex- 
ample, with K = 0.54, 0.71, 0.70, respectively, all of them for laminar flow. The 
experimental constants K come very close to that in our eq. (39a) with k = 1.0, but 
the range of Ra is similar. The details are shown in Fig. 5 • This discussion ap- 
plies to the average Nusselt numbers only. It is seen that eqs. (2l) (horizontal 
plate heated on one side only), and (24) (horizontal plate heated on both sides), 
this paper, fall unquestionably into the range of the experimental results reported 
in the recent literature, although there are differences in the power of the Ray- 
leigh number: m = l/5 is recommended in reference 14 and in this paper (depending 
on the flow situation), and in reference 14, while references 15 to 18 report m = 
l/4. Also, in this paper, expressions for m = l/4 and m = l/3 have been derived. 

For turbulent convection results shown in Fig. 4, the least squares fit equa- 


51 



tion is 


Nu l = 0.0799 Ra^ 0,384 (49) 

that agrees reasonably well with the calculated formula, eq. ( 32 ) here, while in re- 
ference 14 the experimental result 

Nu = 0 .173 Ra 1 ^ 3 (50) 

r g 

is proposed, valid for 10 ^ < Ra < 10 . This result is typical of other experimental 
findings shown in Fig. 5, where, in addition, references 14 to 21 are included. 

From the review of the above it is seen that m = l/3 is very common for the situa- 
tions where turbulent free convection takes place on horizontal plates, heated on 
one side only. For thin plates, of finite size, heated on both sides in finite en- 
closures, m increases somewhat, due apparently to flow instabilities generated by 
heat transfer at the bottom of the plate . 

To produce the exponent m = 1/3 for turbulent, free convective flow in the 
present case, p = 1/3 would have been necessary in eq. (25). It is interesting to 
see that most of the experimental results for the upward-pointing, heated plate 
show a remarkable similarity with the results for the vertical walls; the finding 
of Eckert and Jackson in reference 4 of m = 2/5 is worth noting, as it is derived 
from a semi-empirical theory, as exemplified here by eqs. ( 25 ) and ( 26 ), and is al- 
so backed up by experimental data. It appears that this similarity in free con- 
vective heat transfer, in comparing results for horizontal with those for the ver- 
tical geometry, lies in the effect of the thermal jet generated at the point where 
the boundary flow meets head-on and a thermal plume is formed - something like the 
reverse of the situation depicted in Fig. 3- This aspect has also been analyzed 
above separately for laminar flow. It is mentioned also by Bosworth in reference 1 ?. 


For a horizontal plate facing downward, the controversy concerning the vali- 
dity of the exponents m= l/5 vs. m = l/4 can be resolved with the analysis of the 
flow situation, represented here schematically as Fig. 3* The flow on a downward- 
facing, heated horizontal plate is visualized. It is seen that near the stagnation 
point a thin boundary- layer- like region exists. Then, farther away, the "boundary- 
layer" becomes at first thicker, but that thickness is gradually diminishing, as 
the edge of plate is approached. In this fashion, two kinds of existing analytical 
solutions, and the corresponding existing dual experimental correlations could be 
reconciled, as either of the flow modes may predominate in the given experimental 
set-up. 


Our own experimental data are shown in Fig. 6. It is seen from Fig. 6 that the 
flow regime there must have been already turbulent for Ra as low as 6 x 10- 5 . There 
is a considerable scatter of the data, indicating a highly unstable flow. Heat 
transfer intensity, as represented by the corresponding Nusselt numbers, is roughly 
50 % of that existing on the upward-facing side of plate, which is also to be found 
in the literature (cf. reference 19). The results shown in Fig. 6 may be correlated 
by either of the two equations below. Equation (5l) represents the least squares 
fit 


Nu L = 0 . 001 7 Ra L ° * 596 ( 51 ) 

while eq. (52) represents the limiting case of flow corresponding to p = 0 in eq. 
(25). or 

Nu L = 0.0064 Ra L °* 5 (52) 
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The relatively large exponents of the Rayleigh number in the present case in- 
dicate that some recirculation of the air in the test chamber was actually taking 
place, as the result of simultaneous heating of the bottom and the top of the test 
plate, thus representing a kind of synergistic effect. Few other data for the pre- 
sent configuration are available for the turbulent regime, for the sake of compari- 
son, except those in reference 19 > correlated by a formula that represents simply 
one-half of the results obtained with the heated flat plate facing upward, with m = 
l/ 3 , within ± 10 ^, as 

Nu=0.06Ra 1 / 3 (53) 

Obviously, more data for this configuration are still required, and our objective 
has been to fill this need at least partially. Additional cases for this geometry 
are discussed in reference 3* A general comparison of results of various investi- 
gators is shown in Fig. ?. It is seen from the discussion above that for laminar 
free convection, on both sides of the horizontal plate, many investigators still 
use correlations of the type Nu = K Ra 1 * while some prefer the supposedly more "theo- 
retical" version of Nu = K Ral/-5. A plausible explanation was given for this discre- 
pancy. In reference 1 ?, K occurring in eq. ( 53 ) is given as 0 . 08 , without details. 

It is worth noting that Mangier' s transformation changes relations valid for 
two-dimensional flow over a wedge into expressions valid also in three-dimensional 
(stagnation) flow, without the change of power of Re, for forced convection. A sim- 
ilar phenomenon, by analogy, may be also expected for free convection, with the Rey- 
nolds number replaced by its physical equivalent, the Grashof number to one half 
power. Moreover, the precise values of K would depend on experimental conditions 
and refinements observed while taking measurements, on geometry, boundary conditions 
on edges and the details of how the plates were heated. It should also be kept in 
mind that the left-hand side of the integrated energy equation is the same for both 
forced and free convection, regardless of the flow regime. The right-hand side of 
that equation must be changed to an empirical formula when the flow is turbulent, 
however . 


CONCLUDING REMARKS 


For turbulent, free-convective heat transfer, characterized by an exponent m> 
0.2 5, most investigators propose m = l/3« It is conceivable, however, that for 
cases of really well-developed, rough, non-isotropic turbulence, that exponent may 
get as high as 0.384 and perhaps even higher than 0 . 5 , as is indicated by the pres- 
ent experiments, for the significantly unstable flow at the bottom of a heated, flat 
horizontal plate, heated on both sides. This uncertainty is also due to the fact 
that experimental results for turbulent free convection for that particular geome- 
try are still very scarce, and additional data are needed. Because of relatively 
low convective heat transfer coefficient here, the effect of radiation correction 
becomes very important . Gases have been found in our experimentation where this 
correction for some specimens amounted to more than 50 % of the total heat trans- 
ferred. Because of its emphasis on the radiation aspects in this case, and since 
both sides of plate were heated, the present investigation may be considered the 
first of its kind. The maximum total experimental error is calculated as ± 15 %. 

A general comparison of the various experimental and analytical correlations 
is shown in Figs. 5 and 7, to stress the degree to which the results from the lit- 
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erature still differ from each other. It is also believed that, for the space en- 
vironment, the present experiment shows that radiation cooling of electronics is 
more promising than free-convective cooling, where characteristic velocities are pro- 
portional to one half power of the Grashof number. With weak gravitational accel- 
eration, Grashof numbers will be small, and effects of convective cooling negli- 
gible . 
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Figure 1. Experimental test set-up. 
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Figure 2. Emissivities of all samples tested. 1 - oxydized aluminum 
with Parson’s black coating; 2 - oxydized aluminum with finger 
prints; 3 - oxydized aluminum, clean; 4 - oxydized aluminum, 
burnished; 5 - polished copper (all results by least squares fitting 
method; error margin: less than ± 6.5%). 


57 





Cl = ST AGN AT I ON POINT FLOWjb=B.L. FLOW 
HEATED PLATE 

/ / y // / >" //////t 



Z 

PLUME 
FORMING 


Figure 3. Schematic representation of natural convection near the 
center of a downward-facing, heated plate of finite dimensions. 



Figure 4. Average Nusselt number vs. Raleigh number for 
three specimens on horizontal, heated plate facing upward 
(maximum error ± 15%). 
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Figure 5. Comparison of results of various investigators, for a heated, 
upward-facing horizontal plate. 



Figure 6. Average Nusselt number vs. Raleigh number for 
three specimens on horizontal, heated plate facing downward 
(maximum error: ± 15%). 
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Figure 7. Comparison of results of several investigators, 
for a heated, downward-facing, horizontal plate. 
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